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Abstract 

In a previous paper (J. Phys. A 36, 11807 (2003)), we introduced the 'asymptotic 
iteration method' for solving second-order homogeneous linear differential equations. 
In this paper, we study perturbed problems in quantum mechanics and we use the 
method to find the coefficients in the perturbation series for the eigenvalues and 
eigenfunctions directly, without first solving the unperturbed problem. 
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1 Introduction 



Boundary value problems play an important role in mathematical physics. 
Eigenvalue problems in quantum mechanics represent some of the most im- 
portant mathematical problems in physics. Schrodinger's equation has exact 
solutions for certain special potentials but in most cases approximation meth- 
ods are required. In perturbation theory the Hamiltonian H = Hq + XVp 
has two parts, a soluble term H to which a perturbation XVp is added. The 
wavefunctions and the eigenvalues of the total Hamiltonian are written as a 
perturbation series in the parameter A with coefficients calculated in terms 
of the exact solutions of H . From the practical point of view, the goal of 
this paper is to devise an iterative method to find the coefficients in such 
perturbation series. 

In this work, we will discuss the application [l]-[4] of the asymptotic iteration 
method (AIM) not to solve the problems directly, as we have done earlier, but 
rather to find the perturbation coefficients. Shortly we shall summarize the 
general structure of the method and show how it is applied to boundary value 
problems. The main purpose of the present paper is to calculate the coefficients 
of the perturbation expansion without using the base eigenf unctions of Hq 
explicitly. The differential equation of the original boundary-value problem 
may first have to be transformed into a certain standard form with appropriate 
properties, such as the presence of non-zero terms in the function and in its 
first derivative, or the incorporation of boundary conditions. 

The paper is organized as follows. In section 2, we outline the asymptotic 
iteration method in section 3 we formulate the application of the method 
to perturbation expansions. In section 4 we apply the theory to a variety of 
interesting problems, and we summarize our conclusions in section 5. 



2 The asymptotic iteration method (AIM) 

In this section we give the main structure of the AIM, details can be found 
in [1]. Given that X (x) and s (x) are sufficiently many times continuously 
differentiable, the second-order differential equation 

y" = X (x)y' + s (x)y (1) 
has a general solution 

X X t 

y(x) = exp(- J adt) C 2 + C x j exp ( J (A (r) + 2a(r))drjdt , (2) 



2 



if for some n > 0, 



s. 



n 




(3) 



A 
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where 



A 



n 




s„_i + A A fc _i, and s n = s^_ x + s A„_i. 



(4) 



The termination condition of the method, as given in (3), can also be written 
as follows 



This 5(x) function has a crucial role for boundary-value problems: the co- 
efficients s n and A„ now depend on the (unknown) eigenvalue E, which is 
determined by the vanishing condition (5). If, for a suitable choice of E, this 
equation is satisfied at every x point, then the problem is called 'exactly solv- 
able'. All the applications in the present paper have this feature in the sense 
that we find exact expressions for the perturbation coefficients. The asymp- 
totic iteration method was developed with the idea of using minimal algebraic 
computations to solve second order differential equations. Indeed, with the 
availability of symbolic mathematical software, the computations required by 
the asymptotic iteration (3) present few difficulties, even for the higher it- 
eration steps. AIM has proved to be useful in obtaining solutions for many 
boundary- value problems, including eigenvalue problems of the Schrodinger 



3 The application of AIM to perturbation problems 

3. 1 Perturbation expansions for eigenvalues 

The purpose of this section is to show that the asymptotic iteration method 
can be used to calculate the coefficients of a perturbation expansion for Schrodinger 
eigenproblems. We suppose that we have a potential that is the sum of two 
parts expressed by 



where V\{x) has an exact solution, V^x) is the perturbation, and A is a pertur- 
bation expansion parameter. We suppose that the eigenvalue can be written 




(5) 



type [l]-[3]. 




(6) 
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as a perturbation series of the form 

E = E<® + XE^ + X 2 E^ + X 3 E^ + 



(7) 



Our goal is to calculate the coefficients E^\ where j = 0, 1,2, . . .. For the 
potential given in (6), the correspondence Schrodinger equation has the form 

Id 2 \ 

1 ' + V 1 (x) + XV 2 {x) = E^{x). (8) 



dx 2 



After writing ip(x) = yo(x)f(x) (and, perhaps, using additional changes of 
variable), we obtain a second order linear homogeneous differential equation 
for the factor f(x) with the general form 

f"(x) = X (x,X)f'(x) + s (x,X)f(x). (9) 



We assume that E has the form given in (7). Consequently, the X (x, A) and 
s (x,X) functions depend on each E^>) too. We now apply the asymptotic 
iteration method to (9). After calculating the X n (x, X) and s n (x,X) functions 
by using (4), we can construct the S(x, X) function from (5) as follows 

S(x, A) = s n (x, X)X n+1 (x, A) - s n+1 (x, X)X n (x, A) = 0. (10) 



If we expand 5(x, X) about A = 0, we get the following series 



5(x,X) = 5(x,0) + 



A 35(x,X) 



1! dX 



+ 



X 2 d 2 5{x, X) 



A=0 



2! OX 2 



+ 



A 3 d 3 5{x, A) 



A=0 



3! OX 3 



+ ... 



A=0 



= Y / X k 5 (k \x). 

According to the procedure of AIM, 5(x, A) must be zero; if this is to be so 
for every A value, then every term of the series must be zero. That is to say, 



5 U) (x) - I ^,A) 
° [X) ~j\ dXi 



0, J =0,1,2, 



(12) 



A=0 



It is clear that solution of the equation 6^(x) = gives us E^ and 6^(x) = 
gives E^\ the first correction term to the eigenvalue. This procedure thus 
allows us to find the coefficients of the eigenvalue expansion. It is an attractive 
feature of AIM that we can also calculate the eigenfunctions, as we now show. 
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3.2 Perturbation expansions for the wave function 



When the condition given in (5) is satisfied, the solution of the (9) can be 
found using (2). If we rewrite the first part of (2), we have 

f( x ) = C 2 exp (- fa(t,X)dt) , (13) 



where a(x, A) = f^f 2 ^, k is the iteration number, and C 2 is an integration 
constant. Again if we expand a(x, A) about A = 0, we obtain the following 
series 



a(x, A) = a (0) (x) + Aa (1) (x) + A 2 a (2) (x) + A V 3) (x) + . . . 

oo 

= Y,^ (k) (x), (14) 

k=0 

where a^\x), j = 0, 1, 2, ... is given by 

a U) (x) 1 (* (s k (x,X)\\ _ 1 ( &a(x,\) \ 

The function f(x) in (13) becomes 

oo 

f(x) = C 2 fM(x)fW(x)fW(x)fW(x)... = C 2 J] (16) 

fc=0 

where 

fW(x) = exp ^-A fc J(a (fe) (0)^ , fc = 0, 1, 2, (17) 

4 Some applications 



In this section, we present some examples to demonstrate how AIM can be used 
to determine the perturbation coefficients. Exact results are recovered for some 
well known cases. We shall discuss three examples: the quartic anharmonic 
oscillator, the complex cubic oscillator, and a perturbed Poschl- Teller poten- 
tial. Next we shall obtain a perturbation expansion for the angular spheroidal 
eigenvalues by using the expansion for the perturbed Poschl- Teller potential. 
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For perturbed harmonic oscillators we note that exactly similar treatments 
may be given for the wider class of potentials of the form V(x) = x 2 + Xx q , 
where q = 3,4,5, .. .. We shall first look at two specific regimes of the quartic 
example q — 4. 



4-1 The quartic anharmonic Oscillator: 



Case 1 



The potential for the quartic anharmonic oscillator may be written V(x) = 
x 2 + Ax 4 . Several interesting studies of this problem may be found in the 
literature [5]- [18]. With this potential, Schrodinger's equation becomes 

/ ft 2 \ 

' ' +x 2 + Xx A )tP(x) =Eip{x), (18) 



alx 2 



where A is a small parameter. After substituting ip(x) = exp( — ^x 2 )f(x) into 
(18), we get 

f"{x) = 2xf'(x) + (1 - E + Xx 4 )f(x), (19) 



where we suppose that 

00 

E = Y / ^ k E (k) - (20) 

k=0 



Now we can calculate the coefficients in the energy expansion. If we solve (12) 
for j = together with (19) and (20), we immediately find 

E^ = 2n+1, n = 0,1,2... 

For j — 1 we find that 

(1) _ 3 15 39 75 

~ 4' T' T' T"" 



Consequently, 
For j = 2, we find 

E (2) 



E W = ^2n 2 + 2n+l). 
4 



21 165 615 1575 



16 16 ' 16 ' 16 
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that is to say 

E {2) = -j-{34n 3 + 51n 2 + 59n + 21). 
16 

Calculating in a similar way, we find that 

(3) _ 333 3915 20079 66825 171153 369063 
n ~ ~6T' ~6T' 64 ' 64 ' 64 ' 64 ' 

and 

E {3) = ^-{125n 4 + 250n 3 + A72n 2 + 347n + 111). 
64 

Finally, we can write E n , for n — 0, 1, 2, . . ., as follows 



£ n = (2n + 1) + ^(2n 2 + 2n + 1)A - y^(34n 3 + 51n 2 + 59n + 21)A 2 
+ ^(125n 4 + 250n 3 + 472n 2 + 347n + 111)A 3 + . . . . (21) 



Case 2 



The Rayleigh-Schrodinger perturbation expansion is not applicable for suf- 
ficiently large value of the parameter A. This drawback can overcome using 
AIM by means of a substitution of the form ip(x) = exp( — |ax 2 — \(3x A ) j \x) . 
After inserting this expression for i/j(x) into the (18), we obtain 

f"(x) = 2(ax + (3x 3 )f'(x) + {a-E- (3 2 x 6 )f(x), (22) 



where a 3 — a — |A = and f3 — We point out that when A > 0, then 
a > 1. If we change the independent variable x to z, where x — then we 
find (in terms of z) 

f"(z) = 2(z + V z 3 )f'(z) + (e - V 2 z*)f(z), (23) 



where t] = = ^r, and e = 9L ^-- It is clear that rj is always less than 1. We 
can write e as follows 

oo 

e = V k t (k) (24) 

k=0 



We now apply our perturbation procedure to this problem. If we solve (12) 
for j = 0, together with (23) and (24), we find e( ) = — 2n. For j — 1 we find 
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that e« = -3n 2 . For j = 2, = |(68n 3 + 30n 2 + 64n + 15). Finally, we can 
write E n by using E n = (1 — e n )a and r\ = ^ as follows 

. 3n 2 . 68n 3 + 30n 2 + 64n + 15 x2 . . 

^ = (2n + l )a + — A _ A 2 + .... (25) 



It is clear that, (25) is more strongly convergent than (21) because (25) in- 
cludes the a term which is always bigger than 1, if A > 0. 



4-2 The complex cubic anharmonic oscillator 



In this section, we will derive a perturbation expansion for the complex cubic 
anharmonic oscillator. This example of a PT-symmetric potential with real 
eigenvalues has had much attention in the literature [19]- [29]. Specifically, we 
consider Schrodinger's equation for the potential 

V(x) =x 2 + i\x 3 , (26) 
namely 

("^ + X2 + tXx3 ) = E ^ X) - (27) 

The eigenvalues E are known to be real and positive. If we write the wave 
function ip(x) = exp( — ^x 2 )f(x) and substitute this form into (27), we find 

/" = 2xf + (1 - E + i\x z )f. (28) 
For small A, the eigenvalue E can be written as the perturbation series 

oo 

E = Y / ^ E(k) - (29) 

k=0 



If we solve the (12) for j = 0, together with (29) and (28), we immediately 
find E@> = 2n + 1. For j = 1 we find that Eg* = 0. For j = 2, we find 

E {2) = 30n 2 +30»+ll ^ Calculating £(3) m & similar we find that #(3) = Q 

Meanwhile for E^\ we must solve the equation for j = 4 and we find 

= — f94n 3 + 141n 2 + 109n + 3l) .. 
256 V ) 

Finally, we can write E n as follows 
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E n = (2n+1) + f jA 2 

15 

+ 7^ f 94 ™ 3 + 141 « 2 + 109 ™ + 3l) A 4 + . . . . (30) 
256 v ' 

This expansion is in agreement with [30]- [32] and same as the result of the 
standard Rayleigh-Schrodinger perturbation expansion. 



4-3 Perturbed Poschl- Teller potentials 
Consider Schrodinger's equation 



-tf'(x)+V(x)1>(x)=Eil>(x), V(x) = a(a + ^ +Acos 2 x, (31) 

sin x 



where < x < ir. After substituting ip(x) = sin a+1 (a;)/(a;) and making the 
convenient change of variable y = cosx, we find the following second order 
differential equation for the function f(y) 

fM2a+3) ^ + ^r-E^yy . (32) 

We assume that the eigenvalue E has the expansion 

oo 

£ = ]TA fc £ (fc) . (33) 

fc=0 

If we consider 8^\x) = for this problem we find E^ = (a + n) 2 , where 
n — 1, 2 .... If we solve the equation 8^(x) = 0, we obtain 

E (i) = (2n-l)q+(n-l)(n+l) 
2(a + n-l)(a + n + l) 

To obtain the second correction to the eigenvalue, we have to solve the equation 
£( 2 )(x) = 0, and from this we find 



E^ 



Pn{a) 



16(a + n - 2)(a + n - If (a + n + If (a + n + 2) ' 
where 

Pn( a ) — u 5 a5 + u 4 a4 + U Z°? + u 2& 2 + Ul a + "^0; (34) 
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and for n — 1,2,3, . 



u = -(n-2)(n- lf{n + l) 2 (n + 2) 

Ul = -3(n - l)(n + l)(n 3 - n 2 - 3n + 1) 

W2 = -i(9n 4 - 24n 3 - 19n 2 + 24n + 10) 
2 

u 3 = 2n 3 + 13n 2 - 5ra - 13 
u 4 = 8n 2 + 2n - 9 
w 5 = 2(2n - 1). 

Thus we obtain the perturbation series for the eigenvalues as follows 



^ =(o+n) > + < 2 »- 1 >«+(»- 1 ><»+ 1 ) 



2(a + n- l)(a + n+ 1) 



A 



PJa) 



16(a + n - 2)(a + n - If (a + n + l) 3 (a + n + 2) 



A 2 + 



(35) 



^.^ TTie angular spheroidal eigenvalues 



From (35), we can get an explicit formula for the angular spheroidal eigen- 
values. This important problem has a long history [33] -[42]. The angular 
spheroidal wave equation can be written in following form [4] 



/" = (2m + 2) -H-^f + ^— - f, (36) 




with — 1 < y < 1. Here, e = (Aj"(c)) 2 — m(m + 1), / is the angular momentum 
quantum number, and m is the eigenvalue of the operator L z . We note that 
for c = 0, (A™(c)) 2 = 1(1 + 1). When we compare (36) with (32), we find that 
a = m — |, A = c 2 , and e = £"„ — (a + l) 2 . After substituting n = I — m + 1 
and a = m — | into the (35), we obtain the perturbation expansion for A™(c) 
shown below: this agrees with the results of [40] for the special case m — 0. 



(Arw) 2 =i(i+i) + 



(2/ - 2m + l)(2m - 1) + 2(Z - m)(l - m + 2) 



Pi(m) 



(21 - l)(2Z + 3) 



2(2/ -3)(2/- 1) 3 (2/ + 3) 3 (2/ + 5) 



e + . . . 



(37) 



10 



where 

P^m) = W5 (2m-l) 5 +2M4(2m-l) 4 +4M 3 (2m-l) 3 +8M 2 (2m-l) 2 +16Mi(2m-l)+32wo, 
and 

u = ^(t-l)t 2 (t + 2) 2 (t + 3) 

Ul = -3t(t + 2)(t 3 + 2t 2 -2t-2) 

u 2 = — ^(9t 3 + 12t 2 - 37t - 50) 

u 3 = 2t 3 + 19t 2 + 27t-3 
u 4 = 8t 2 + 18* + 1 
u 5 = 2(2t + 1) 

where t — l — m. 



5 Conclusion 

The asymptotic iteration method allows one to iterate a second-order linear 
differential equation so that at each iteration an equation of the same general 
type is recovered. This iteration scheme applies to boundary value problems 
provided the boundary conditions are accommodated by a factored form for 
the wave function that is asymptotically correct. An exact eigenvalue is ob- 
tained when the energy E can be chosen so that the iteration converges in a 
finite number of steps. In more general cases an approximation is obtained by 
forced numerical convergence. In this paper we have applied the method to 
some perturbation problems and have been able to find the coefficients in the 
perturbation series directly, without first solving the unperturbed problem. 
In many cases, such as those studied here, we have been able to derive these 
coefficients exactly. 
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